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Motivation 


* Monotonic Basin Hopping (MBH) has been shown to be very effective 
optimization methodology for interplanetary space-probe trajectories — which are 
often very challenging constrained global optimization problems in high- 
dimensional solution spaces. 


* When trajectory optimization methods are autonomous, “when to stop” 
optimizing is an important question. 

Despite the many benefits of MBH, we do not see how it can autonomously 
determine when the global search is complete. 


* The goal of Filament Walking (FW) optimization is to address the “when to stop,” 
which is related to “thoroughness” in finding and walking filaments. 


* This presentation is a progress report. It does not present the solution to the 
“thoroughness” issue, or — thereby — the “when to stop” question. 


It presents FW as an alternative to MBH, demonstrates that the assumptions of 
FW are realistic and it’s logic is sound, brings new precision to the “when to stop” 
question, and outlines future work. 


The basic concept of FW 


Optimization of interplanetary space-probe trajectories is a difficult and complicated 
constrained global optimization problem in high-dimensional solution spaces. 


Such problems are more challenging when low-thrust propulsion is used. 


Often, the massive optimizations that are required must be totally autonomous. They do 
not allow time for human intervention (e.g., for initial guesses). 


We “shrink” the high-dimensional problem space and its complicated constraints, down 
to a small number of 1-dimensional search with simpler constraints. 


We do this by transforming the set of candidate-solution points, that are feasible with 
respect to equality constraints, into one (or a few) one-dimensional search filaments that 
can be found and walked along “thoroughly” with extremely high probability. 


At each step in the walk along a filament the objective function at that point is sampled. 
Then we applying the equality constraints to the points on the filament. 


Formal problem statement 
* We define x as a point in a dense set of points X in solution space S in RN 
* In typical problems, N is approximately 10 to 100 
* We seek global min(f(x)), subject to: 
= g,(x)-c, = 0 for all | equality constraints 
= h,(x)-k; < 0 for all J inequality constraints 
= In other words, we seek global min(f(x**)), where x** are strictly feasible points 


= We say that points x* are feasible at least with respect to the set of all equality 
constraints g,(x)-c; 

= Because we require the x* to satisfy all of the equality constraints, they can only reside at 
the intersections of the g,(x)-c, = 0. We call these intersections “filaments.” 


= Because the functions g,(x) that determine the equality constraints are known to be first- 
differentiable, and they are usually are reasonably “smooth”, the gradients of g;(x) near 
an x* usually “point” to that x*. 

" The x** may then be found by filtering the x* to remove all cases where h,(x*)-k; > 0. 


= x* is then found by evaluating all of the f(x**) 


Our initial concerns about filaments 
Their existence, their lengths, and the number of them 
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The apparent nature of real problems 


* So far, filaments do not seem to be too “short” or numerous in real problems. 
* Gradients do not seem to be “rippled,” “rugged,” roughly textured, or locally flat. 
* We are left with the challenge of “thoroughness:” 

* How can we know that we have “almost surely” found all the filaments? 


* How can we know that the walker has walked, “almost surely,” each filament 
thoroughly enough to have found that filament’s min(f(x*))? 


* Can we build autonomous quantitative measures of “thoroughness?” 


Simplified example 


Wireframe of objective function f 


Heatmap of f with filament shown 
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Testing our simple example 


5,000 steps along the filament 


min(f(x*)) has been found (twice in 5,000 steps) 


f along the filament 


+} Length of the filament —— 


True (a prior known) min(f) = -2646.099011 
True (a prior known) argmin_x(f) = 500 on x = (0, 1000] 
True (a prior known) argmin_y(f) = 500 on y = (0, 1000] 


f along the 5,000-step filament walk 
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min(f) along the filament found by the walker = -2646.099011 
argmin(s)_t(f along the filament) = 1742 and 3251 on t = (0, 5000] 
argmin_x(f along the filament) = 500 on x = (0, 1000] 

argmin_y(f along the filament) = 500 on y = (0, 1000] 


}<—_——Time (steps) spent walking 
the filament 


Success with the simple example 
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After 5,000 steps along the filament 
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Applicability: Our initial concerns and what we found 


* Do filaments really exist in important trajectory optimization problems? Yes 
* Can they be found? Yes 

* Are they walkable? Yes 

* Are they few in number? So far, there at most a few 


* Are they unique? Not always 


Are they sufficiently “long,” so as to be “walk-worthy” Yes, so far 


* Are the gradients that “point” towards them sufficiently smooth that if the walker gets 
“close” to a filament, or falls not-too-far off of it, the gradients will “guide” it back to the 
filament? Apparently yes. (Much better than one of us thought) 


If filaments exist, can be found, are long, and are walkable, how will we visualize them? 
You will see in the next slides 


Context for the real problem we are about to show: 
Multiple Shooting Transcription for High-Thrust Chemical Propulsion 


| Mateh point @ Control point (planet, satellite, free point in space) 


f Deep space maneuver 


Forward propagation 


Mission clock time flow 
—_—> 


m,=m, 


Courtesy of Donald Ellison and Matthew Vavrina 


Earth to Mars Transfer Trajectory Example 


Transcription MGAnDSMs* [10] 
Launch date bounds 2/16/2006 - 11/12/2008 
Departure type launch with C, < 48.5 km2/s? 
Arrival type chemical rendezvous 
Objective function minimize Av 


Flight time upper bound 500 days 


Known global optimum 2.343 km/s 
(from MBH) 


*Multiple Gravity Assist with n Deep Space Maneuvers using shooting 


Image of a filament from the first real problem 


* Aprojection of a 13-dim 
space onto to two 
dimensions (normalized 
launch date versus 
normalized flight time) 


* This is the filament walking 
equivalent of the so-called 
“pork-chop” plot 


normalized flight time 


Objective function 


normalized launch date 


Unraveling the filament 


A distance-preserving homeomorphism of the walk 


The above suggests: 


Only one filament exists 


It was walked reasonably thoroughly, mostly in a single direction 
Occasionally, the walker “turned around” (switched directions) 
When the walker fell off a filament, it was usually guided back, quickly, onto that (or 


another) filament by local gradients 


Homeomorphism 


Future work 


Can we assure that all of the filaments were be found “almost surely?” 


Can we assure that a walker has walked each found filament “thoroughly” 
enough to find min(f(x*)) on that filament “almost surely?” 


Can we predict the number and “length” of filaments that should be to expected, 
based on the structure of the problem? 


Can we build a autonomous quantitative measures of confidence in whether the 
true global min(f(x*)) has been found? 


Can one design transcriptions that are especially well-suited to filament walking? 


Benchmarking FW against MBH in terms of accuracy, speed, and reliability 


Summary 


* The concept and assumptions of FW appear to be valid and applicable. 
* The “when to stop problem” is not yet solved but, for the first time, we have a 
promising framework in which to solve it. 


* This is the first phase of a long journey. 


Thank you 


Arnold C. Englander 
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